arXiv:1509.05847v3 [cond-mat.mes-hall] 12 Apr 2016 


Spin torque and Nernst effects in Dzyaloshinskii-Moriya ferromagnets 
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We predict that a temperature gradient can induce a magnon-mediated intrinsic torque in systems 
with non-trivial magnon Berry curvature. With the help of a microscopic linear response theory of 
nonequilibrium magnon-mediated torques and spin currents we identify the interband and intraband 
components that manifest in ferromagnets with Dzyaloshinskii-Moriya interactions and magnetic 
textures. To illustrate and assess the importance of such effects, we apply the linear response theory 
to the magnon-mediated spin Nernst and torque responses in a kagome lattice ferromagnet. 
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Studies of the spin degree of freedom in spintronics 
[Tj naturally extend to include the interplay between the 
energy and spin flows in the field of spincaloritronics 
H E]. Improved efficiency in interconversion between 
energy and spin [4] could result in important applica¬ 
tions, e.g., for energy harvesing, cooling, and magnetiza¬ 
tion control isHini. Magnetic insulators such as yttrium 
iron garnet (YIG) or LU 2 V 2 O 7 offer a perfect playground 
for spincaloritronics where due to the absence of electron 
continuum the dissipation can be lowered as only the spin 
and energy matter HIHISI- It has already been demon¬ 
strated in recent experiments that energy currents can be 
used for magnetization control [HUS]. This opens new 
possibilities for applications of magnon-mediated torques 
in racetrack memories dsiin], and even in quantum in¬ 
formation manipulations |18| . 

As we show in this study, the magnon-mediated torque 
is closely related to the magnon-mediated thermal Hall 
effect. The latter has been observed in LU 2 V 2 O 7 m 
and explained by the Berry curvature of magnon bands 
usHm where the physics is reminiscent of the anoma¬ 
lous Hall effect [H]. The possibility of the magnon edge 
currents and tunable topology of the magnon bands has 
also been discussed in the context of magnetic insula¬ 
tors mi In a recent experiment, the magnon- 

mediated thermal Hall effect showed the sign reversal 
with changes in temperature or magnetic field in the 
kagome magnet Cu(l-3, bdc) [26]. Since magnons also 
carry spin it would be natural to also study how spin 
currents can be generated from temperature gradients, 
i.e., the spin Nernst effect, in materials with nontrivial 
topology of magnon bands. However, both the magnon- 
mediated torque and the spin Nernst effect have not been 
addressed in systems with non-trivial magnon Berry cur¬ 
vature. Such calculations inevitably require generaliza¬ 
tions of linear response methods developed in sixties and 
seventies EH HE] to bosonic systems and consideration 
of the spin current analog of the energy magnetization 
contribution |29| . 

In this Rapid Communication, we predict that a tem¬ 
perature gradient can induce a magnon-mediated intrin¬ 



Figure 1. (Color online) Left: Two-dimensional kagome lat¬ 
tice with lattice vectors ai and a 2 where atoms are placed 
in the corners of triangles. DMI are shown by vectors per¬ 
pendicular to the bonds. The DMI of strength Di point into 
the page while the Rashba-like DMI of strength D 2 lie in the 
page. Right: The three magnon bands are plotted for the 
case of broken mirror symmetry with respect to the kagome 
plane due to the Rashba-like DMI. The direction of the spin 
density is given by m = :rsin(7r/6) + zcos(7r/6). 


sic torque in systems with non-trivial magnon Berry cur¬ 
vature. To this end, we formulate a microscopic lin¬ 
ear response theory to temperature gradients for ferro¬ 
magnets with multiple magnon bands. We follow the 
Luttinger approach of the gravitational scalar potential 
[271130]. Our theory is capable of capturing the nontrivial 
topology of magnon bands resulting from the Dzyaloshin¬ 
skii-Moriya interactions (DMI) [311 [32]. An additional 
vector potential corresponding to the magnetic texture 
can be readily introduced in our approach via minimal 
coupling. We note that the predicted magnon-mediated 
torques are bosonic analogs of the spin-orbit torques [33- 
|42] . We find that torques due to Dzyaloshinskii-Moriya 
interactions (DM torques) can only arise in systems lack¬ 
ing the center of inversion. This is in contrast to the 
the magnon-mediated spin Nernst effect. Finally, we ex¬ 
press the intrinsic contribution to the DM torque via 
the mixed Berry curvature calculated with respect to the 
variation of the magnetization and momentum [22]. We 
apply our linear response theory to the magnon-mediated 
spin Nernst and torque responses in a kagome lattice fer- 
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romagnet. We note that the latter can be detected by 
studying the magnetization dynamics while the former 
can be detected by the inverse spin Hall effect. 

Microscopic theory. — We consider a noninteracting 
boson Hamiltonian describing the magnon fields: 

H = J dr^\r)H^{r), (1) 

where is a Hermitian matrix of the size N x N and 
= [a|(r),..., a]y(r)]) describes N bosonic fields 
corresponding to the number of modes within a unit 
cell (or the number of spin-wave bands). Hamiltonian 
in Eq. 0 can also account for smooth magnetic textures 
via minimal coupling to the texture-induced vector po¬ 
tential A via additional term {A.a • where is the 
magnon spin current [Mills]- The Fourier transformed 
Hamiltonian is: 

n = J2'^iHi^)ak, ( 2 ) 

k 

where a\ is the Fourier transformed vector of creation 
operators. Hamiltonian in Eq. 0 can be diagonal¬ 
ized by a unitary matrix T/^, i.e. = TlH{k)T]^ and 
T^T/c = InxN where Sk is the diagonal matrix of band 
energies, and InxN is the N x N unit matrix. As it was 
shown by Luttinger mi, the effect of the temperature 
gradient can be replicated by introducing a perturbation 
to Hamiltonian in Eq. 0 : 

'H' = \ I dr^\r) {Hx + xH) ^(r), (3) 

where the nonequilibrium magnon-mediated field can be 
treated as a linear response to the perturbation in Eq. ^ 
and diX = diT/T. 

The nonequilibrium magnon-mediated field can be cal¬ 
culated by invoking arguments similar to those for the 
spin-orbit torque [sasnss]: 

htot = h + h = (4) 

where the averaging is done either over the equilibrium 
or nonequilibrium state induced by the temperature gra¬ 
dient, and m is a unit vector in the direction of the 
spin density s. The magnon-mediated torque can be ex¬ 
pressed as T = m X htot leading to modification of the 
Landau-Lifshitz-Gilbert equation, i.e., 5(1 + (amx)m = 
m X Heff + r where Heff is the effective magnetic field. 
We are also concerned with the magnon current carrying 
spin which has two components: 

Jtot = (JL + (j'),,, (5) 

where the first component, J = f (irT’t'(r)vT(r), does 
not depend on the temperature gradient and the second 


component, j' = (1/2) J drT’t'(r) (vy + yv) T(r), is pro¬ 
portional to the temperature gradient. The latter contri¬ 
bution is related to the spin current analog of the energy 
magnetization [29]. Here the velocity operator is given 
by V = (l/ih)[Y^H]. The magnon current density j is in¬ 
troduced in a standard way from the continuity equation 
p +V •j(r) = 0 where p is the density of magnons. In our 
discussion, we employ the expression for the energy cur¬ 
rent density, j^(r) = (l/2)T’t'(r)(vi7-f i7v)T(r), and the 
macroscopic energy current = J drj^(r) correspond¬ 
ing to the continuity equation p£; + V- j^(r) =0 with 
pE being the energy density. Note that we omitted the 
component of proportional to diX as it is irrelevant 
to our discussion. Within the linear response theory, the 
response of an operator X to temperature gradient be¬ 
comes: 

A)ne = “ n^(0)]A^} d,x, (6) 

where X is either spin current — fiJ or nonequilibrium 
field h = —dmH^ n^(0) = + iO) is the retarded 

correlation function related to the following correlator in 
Matsubara formalism, Ilij{ift) = — (T^XiJ^). 

Note that the energy current originates from the expres¬ 
sion v! = {ilh)\H, V!] = 

We calculate the correlator in Eq. 0 by considering 
the simplest bubble diagram for H^j and performing the 
analytic continuation. We express the result through a 
response tensor tij = t\j -h such that Xi = —UjdjX 

m- 

= 11 - XiG^JjG^), 

,■■ 4 /^»(«-)EeTV(X.G V, 

(7) 

where g{uj) is the Bose distribution function g{uj) = 
l/[exp{huj/kBT) - 1 ], = h{huj ir)-\ = 

h{huj — H — ir)~^, and J' = {vH + H'v)/2. For practical 
purposes, we Fourier transform Eq. 0 which leads to 
additional momentum integration and momentum trans¬ 
formed terms, i.e. G^(k) = h{huj — H'(k) + 'i^T)”^, 
h/c = v/e = and Jk = (v/,i^(k) + 

H’(k)v/c)/2. The approximation we are using can be im¬ 
proved by performing the disorder averaging which is in¬ 
dicated by brackets in Eq. 0 . In addition, interactions 
with phonons can be also taken into account and can re¬ 
sult in additional dissipative corrections to the torque. 
Throughout this paper, we adopt a simple phenomeno¬ 
logical treatment by relating the quasiparticle broadening 
to the Gilbert damping, i.e. T = ahw. 

Berry curvature formulation. — It is very insightful to 
carry out the frequency integrations in Eq. 0 , keep¬ 
ing only the two leading orders in T and combining the 
linear response result with the nonmuilibrium contri¬ 
bution h in Eq. Q or J in Eq. To carry the 
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integrations in Eq. 0 we use the diagonal basis de¬ 
fined by rotation matrices T/., and transform the con¬ 
tributions h and J to an integral over energies follow¬ 
ing the approach of Smrcka and Streda [iniiiH]. Using 
the covariant derivative we calculate the rotated veloc¬ 
ity, TlhvkTk = d\^Ek - iAk£k + iEkAk, and nonequilib- 
rium field, - iAm^k + i^k^, where 

A-k = i'^dkTk and Am = iT^OmTh- Substituting these 
in Eq. Q we identify the intraband and interband con¬ 
tributions to the response tensor |46| : 

1 ^ 1 

2p {^Xi^nk){^kj^nk)^nk9 {,^nk\ 
k n=l^ (o\ 

h T ^ ^ 

k n=l 

where Xi is either or ki, Snk = [£k]nn, '^nk = a£nk, 
g'i^nk) = {2kBT)~^{l - cosh(e„fc/fcBT)}“\ ci[x] = 
fg dt ln[(l + t)/t] = (1 +a:;) ln[l +a:;] — x Ina:;, V is volume, 
and we introduced the Berry curvature of n-th band: 

1^" (k) = n(d,,Tl)(dk,n) - (dk,Ti)(d,,n)U. (9) 

Such Berry curvatures naturally appear in discussions of 
semiclassical equations of motion for Hamiltonians with 
slowly varying parameters [22] . Derivation of Eq. ^ (see 
supplemental material [46|) should also hold for fermion 
systems given that ci(snk) = ~ dn^(v) where 
the Eermi-Dirac distribution np replaces g m By ap¬ 
plying the time reversal transformation, i.e. k ^ — k, 
m ^ —m, ^ recover 

the transformation properties of tand under the 
magnetization reversal. In particular, it is clear that t 
is even under the magnetization reversal and is divergent 
as r ^0. On the other hand, is odd under the mag¬ 
netization reversal and corresponds to the intrinsic con¬ 
tribution independent of E. In terms of spin torques, the 
former corresponds to the field-like torque and the latter 
to the anti-damping (or dissipative) intrinsic torque. 

Model. — We apply our theory to the magnon current 
and torque response of a kagome lattice ferromagnet with 
DMI (see Eig. [^. The exchange and DMI terms in the 
Hamiltonian are given by mi ED; 

= Si + ^ El^b ■ (Si X S,), (10) 

where J corresponds to the nearest neighbor interaction, 
is the DMI vector between sites i and j (J^ij = 
—Dji). We take the DMI vector to be = Diz for 
the ordering of sites shown by the arrow inside trian¬ 
gles in Eig.[^ Such configuration corresponds to systems 

with the center of inversion. In some cases, we also add a 

Ti 

Rashba-like inplane contribution, = D 2 {zx ij)^ that 
breaks the mirror symmetry with respect to the kagome 
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Figure 2. (Color online) Left: The spin Nernst conductivity 
afj versus temperature T for DMI D 2 = 0 and Di = J/2, J/3, 
and J/6. Right: The thermal torkance corresponding to 
the anti-damping part of the torque versus temperature T 
for DMI Di = D 2 = J/2, J/3, and J/6. Note that the 
temperature range is not limited by the Curie temperature in 
order to show the assymtotic behavior. In both figures the 
direction of the spin density is given by m = £. 


plane where ij is a unit vector connecting sites i and 
j (Dij is shown by arrows in Eig. [^. We also add the 
Zeeman term due to an external magnetic field that fixes 
the direction of the magnetization direction along the 
field. After applying the Holstein-Primakoff transforma¬ 
tion, we arrive at a noninteracting Hamiltonian compat¬ 
ible with Eq. 0- A typical magnon spectrum is shown 
in Eig. [^where the lower, middle, and upper bands have 
the Chern numbers 1,0, and —1, respectively. 


Spin Nernst effeet. — The thermal Hall effect manifests 
itself in the transverse temperature gradient HlEolES]. 
Here we calculate the transverse spin current which can 
be detected, e.g., via the inverse spin Hall effect in a Pt 
contact attached to the sample [48]. The spin Nernst 
conductivity relates the temperature gradient to the 
spin current density, i.e. jf = —hji = —afjdjT where 
each magnon carries the angular momentum —h. Erom 
Eq. (j^ we obtain afj = tij/T with only the interband 
part contributing to a\y Eor a model calculation, we 
consider Eq. (10). The spin Nernst effect can take place 
in systems with the center of inversion, thus the Rashba- 
like DMI described by parameter can be zero. By 
integrating the Berry curvature over the Brillouin zone, 
we arrive at the result in Eig.[2|where al^ is dominated by 
the lowest band in Eig. with the positive Chern num¬ 
ber. Eor a three-dimensional system containing weakly 
interacting kagome layers, we can write = al-/c 
where c oc a is the interlayer distance and a is the lattice 
constant. Given results in Eig. it seems to be pos¬ 
sible to generate a transverse spin current of the order 
of 10“^^J/m^ from a temperature gradient of 20K/mm 
m in three dimensional systems. Spin currents of such 
magnitude are typical for spin pumping experiments [4]. 


Nonequilibrium torques .— To present our results we 
introduce the thermal torkance /3ij that relates the mag¬ 
netization torque to the temperature gradient, i.e. Ti = 
—PijdjT or Pij = miSikitkj/T in terms of Eq. ^ where 
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Figure 3. (Color online) The nonequilibrium magnon- 
mediated torque 'T is plotted on a unit sphere representing 
the direction of uniform spin density m. The temperature is 
T — 2J and the gradient is applied along the —x direction. 
The field-like torque component that is odd in the mag¬ 
netization is plotted on the left and the anti-damping com¬ 
ponent that is even in the magnetization is plotted on 
the right. The field-like component is rescaled by the Gilbert 
damping to match in scale the anti-damping component, i.e., 


eiki is the antisymmetric tensor. We further separate the 
torkance jSij into the field-like part that is odd in the 
magnetization and the anti-damping part Pfj that is even 
in the magnetization. 


To uncover the effect of Berry curvature, we apply 
our theory to the model in Eq. (10). Within our the¬ 
ory the anti-damping component of the torque entirely 
comes from the Berry curvature contribution in Eq. 

The largest component of /3fj corresponding to the tem¬ 
perature gradient along the x—axis, the torque along the 
axis, and the spin density along the 2 ;—axis is plotted 
in Eig. 1^ The temperature dependence of resem¬ 
bles the temperature dependence of the spin Nernst con¬ 
ductivity where we observe larger effect at higher tem¬ 
peratures. Eor a three-dimensional system containing 
weakly interacting kagome layers, we obtain = /3ij/c 
where c is the interlayer distance. In Eig. we plot the 
nonequilibrium magnon-mediated torque separated into 
the field-like and anti-damping parts, 'T = ^ on a 

unit sphere representing the spin density vector m. The 
torque in Eig. can be obtained from phenomenological 
expressions obtained for films with structural asymme¬ 
try along the z— axis nnmi, t{ (X (m X T>i)diT and 
'T^ oc m X (m X D^)9^T, by a deformation not involving 
the change in topology where x and i is either 

X or y. 


A ballpark estimate of the strength of the nonequilib¬ 
rium magnon-mediated torque can be done by consider¬ 
ing only the lowest band in the quadratic approximation, 
i.e., we have H{k.) = hA[ka^m-{Y>a/A —yAa)]‘^/s where 
A is the exchange stiffness, A.a is the texture-induced 
vector potential, s is the spin density, and a tensor = 
Dq, • e /3 describes DMI. After substituting this spectrum 


in the first Eq. we obtain the longitudinal spin cur¬ 
rent = —h} = kB^T[^/7^({3/2)]/{87^‘^Xa) where ( is 
the Riemann zeta function and A = yjKA/sksT is the 
thermal magnon wavelength. The same Eq. @ results 
in the expression for the nonequilibrium field-like torque 
density: 


r^ = [mx(D„M-^„)]j*, (11) 

which agrees with the earlier results obtained for a single¬ 
band ferromagnet [lOlllSHsl]. Here the torque is gener¬ 
ated within the whole volume. This is contrary to the 
conventional spin-transfer torque which is generated only 
close to the interface m- The typical charge current den¬ 
sity = lO^^A/m^ sufficient for the spin-transfer torque 
switching should be compared to 2ejdDlA ^ lO^A/m^ 
where e is the electron charge, D is the strength of 
DMI and d is the width of the magnet. Eor the esti¬ 
mate of the field-like torque, we assume that d = A/D, 
diT = 20 K/mm, and a = 10 “^ ^TE\ . 

Conclusions. — We developed a linear response theory 
to temperature gradients for magnetization torques (DM 
torques). We identify the intrinsic part of the DM torque 
and express it through the Berry curvature. We note that 
similar expressions also arise for the magnon-mediated 
spin Nernst effect. According to our estimates, the spin 
Nernst effect leads to substantial spin currents that could 
be measured, e.g., by the inverse spin Hall techniques [48] 
in such materials as pyrochlore crystals (e.g., LU 2 V 2 O 7 ) 
and the kagome ferromagnets [M1I54] [e.g., Cu(I-3, bdc)]. 
In particular, a voltage should arise in the neighboring 
heavy metal due to the inverse spin Hall effect in full anal¬ 
ogy to measurements of the spin Seebeck effect and spin 
pumping [3|. We also find that the DM torques should in¬ 
fluence the magnetization dynamics in ferromagnets with 
DMI; however, larger temperature gradients (compared 
to 20 K/mm used in estimates [T5|) are required, e.g., for 
magnetization switching Eor the validity of the lin¬ 
ear response approximation the temperature should not 
change much over the magnon mean free path. The DM 
torque can only arise in materials with structural asym¬ 
metry or lacking the center of inversion. Of relevance 
could be jarosites [56] or ferromagnets and ferrimagnets 
containing buckled kagome layers ISH EH]. Our theory 
can be readily generalized to antiferromagnets and fer¬ 
rimagnets, extending the range of materials suitable for 
observation of DM torques. In particular, antiferromag- 
net does not have to have the center of inversion in order 
to exhibit the DM torque provided each sublattice indi¬ 
vidually lacks the center of inversion. 
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A. The Kubo formula for linear response 


We study a response of an operator X to the perturbation described by a Hamiltonian PL'. The total Hamiltonian 
of the system is: 


Utot^n + W, 


(1) 


where Pi = f and PL' = f We assume that the Hamiltonian H does not contain 

particle non-conserving (pairing) terms. The perturbing Hamiltonian could he H' = ^ (xiJ + Hx) which corresponds 
to the temperature gradient diX = diTjT, and H' = ecj) which corresponds to the electric field Ei = —dif) with e 
being the charge. 

We assume that the operator X can in general depend on the gradients. We are going to consider only the responses 
of the global operators A = y f dr^^ (r)X^(r) where we take V = 1 and reintroduce V in the final expressions. Within 
the linear response, we only need to consider the two lowest orders with respect to gradients: 


Let us now write a general expression for the response of an operator A to the perturbation PL': 

(a} = (a<-°'>) , 

\ / ne \ I eq 


( 2 ) 


(3) 


where < .. >ne is the average over the non-equilibrium states due to the gradient terms, and < .. >eq is the average 
over the equilibrium state. The first average can be found from the Kubo formula of the linear response [1]. It is 
given by the correlation function: 




(4) 


where A{t) = and PL'{t) = Y,mk are the eigen-basis representa¬ 
tions, with the field operators satisfying the relation HAtn{r) = We use the relation PL = {i/h)[PL,PL ] = 

J^dx when integrating by parts to obtain: 




dn'{-t). 

dt 


ZCJ - 0+ h 


i r 

Jo 




dt 


iuj - 0+ 


1 


:dx [Ko(c^)-Ko(0)], 


iuj - 0+ 

where we assumed that the perturbation vanishes as t ^ —oo. We also introduces a new correlator: 

Ko(c^) = ([^(i), J^^]) . 

To proceed further in simplifying the correlator, we use the following identities: 

n) — ^mq^np'kl' {Em) [1 + n(K)] - ^mq^np n{E„) [1 + n{Em)] , 


(5) 

( 6 ) 

(7) 

( 8 ) 


where (...) stands for taking an average over the unperturbed state, n{E) = l/(e^^ — 1) is the Bose distribution 
function, and /? = 1/T. After some transformations, we obtain: 


K r. .'I _ ^ n{Em) - n{En) q 

n,m 


(9) 



2 


which is equivalent to the following expression: 

Ko(c.) = I 


nOO 

/ ft'R 

5{E - Em) 

1 

j LLJ-J iLyJZ/J 

J —oo 

_En - E - UO - i0+ 

E - Em - uo - i0+ J 




ran'' nm ’ 


( 10 ) 


The static limit can be calculated by expanding the last expression in uo. The zeroth order term in uo cancels with the 
K(0) contribution to the correlation function which results in the expression: 


Ko(a;) - Ko(0) « I ^ n{E) 


6{E - Em) 


6{E - E^) 


{E-En + iQ+f {E-Em- *0+)" 


/i(o) iQ 

ran nm ’ 


( 11 ) 


This expression coressponds to the one derived by Crepieux and Bruno in [2]. Identities 5{E — En) = 

— (27ri)“^ — G^{E)) and {E — E^Ei^)^) ^ ^re of use in order to express the correlator in 

terms of Green’s functions: 


, 1 __/*+oo 

=—y/ dEn{E) 


(^) I • (12) 


This expression is consistent with the expression (31) in Smrcka and Streda in [3]. After some transformations, we 
recover expression corresponding to Eq. (7) from the main text: 


= /r S A { } - n{E)Re | 


dGl ^ 

dE ^ ^ dE 


AMm]) dx. 
(13) 


Here, we assumed that the boundary terms vanish, and possible singularity at E = 0 can be regularized. 

It is instructive to show that the same expressions are obtained when calculating the response correlation function 
in the Matsubara frequency. The correlator becomes: 




(14) 


where P = 1/T, uj = 27rTm is the boson Matsubara frequency with integer m, and r = it is the imaginary time. We 
then perform integration by parts: 


— / dry- (TrA^°H0)n'{-T)) = ^- (TrA^°H0)n'{-T)) 

Jo dT \ / 'l^m ' ' 


Jo 


/o ’^^m \ dr 

We will use the relation idrU' = {i/h)\H^'H ] = Observing that = 1, we obtain the correlator: 


,^(°)(0)H'(-r)) 1^ = f^T— /t,A(°)(0) 
' Jo \ 


, J_K,o)ax. 

dr / lUm 


Therefore, for ^ 0 we can write 


(A(o.))„, = — (K(0) - K{u;m)) dx = i 

tuJm, 


du>, 


-K{L 0 m)dx 


(15) 

(16) 

(17) 

(18) 


UJm^O 


B. Total current response to temperature gradient 

Magnon current corresponding to p(r) = ^'^(r)^(r) density, is given 

j = §t(r)v§(r), 


( 19 ) 
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where = (l + ^) ^(r), v = ^[H,r] is the velocity and we defined = ^'^(r)v^(r) and = 

l^t(r) (vr^ + r^v) ^(r)V^X. We will be working with macroscopic currents, defined as J = / drj. When calcu¬ 
lating the linear response, we define 

{Ja) = {j}°^) +(jW) , (20) 

where the first term is an average of Ja^ over the non-equilibrium state, which is just a Kubo linear response expression. 
We then write 

= - lim = 5„W/3X, ( 21 ) 

co>—)-0 OiO 

where Ifc^^ = (r) J^(0)^, where cj is a boson Matsubara frequency, T is temperature, and = 

^ / (ir^'^(r) [Hv^ + v^H] ^(r)V^X is the current occuring when deriving the Kubo formula. The other term in the 
definition of the response {J^) has a straightforward definition 

')eq " ^ ^ ^a/3V/3X. ( 22 ) 

We then write 


{Jo) — [^aP T ^ap] ^pX — pX.. 


(23) 


Let us deinfe a unitary matrix that diagonalizes the Hamiltonian. In Fourier space, we write 

T+ifkTk = ^k, 


where Tj^Ik = 1. In the diagonal basis the velocity becomes: 

l^k,a = ^a^k ^•^k,a^k “t“ ^^k*^k,a5 


(24) 


(25) 


where .4ka = iT^9aT\^, and Vk.a is a Fourier transform of the velocity operator Va- 

Let us calculate the Sap- We denote energy of band n as {Sk)nn = ^nk- After some calculations the expression for 
Sap becomes: 


SaP — ^ ^ (f^k,a)nm (^k'^k,;0 T '^k,;0^k)^ 


9 (^nk) 9 (^mk) 


kn 


(^nk ^mk T 

= -^ E(-4 k.,a)nrn{<-^\<.,p)rnn^nk.9 (^nk) T ^ ^ (A-k,a)nm^mk(A-k,;g) mn g (£nk) - {a-H-13) 

L kn kn 

T ^ ^ (^g^nk) (^;g^nk)^nk 9 (^nk) 

kn 

dg5 [77 - £„k] giv) iv + -ffk) (^/jTk)] - (a <-> ^) + ^{dae„k){d 0 enk)e„k^^g'{enk). 

(29) 


(26) 

(27) 

(28) 


kn 


We define two distinct by nature terms: 



5] r dvS iv - £„k) giv) [ (daTl) iv + Hk) idpTk) 

kn 


{a ^ (3), 


and 



^ ^ (^g^nk)(^;0^nk)^nk 

kn 


1 

2r^k 


9 (^nk)- 


(30) 


(31) 


[21 

In deriving the S^^p term we assumed a lifetime of the bosons F^k, and used the transformations from the previous 
section for Eq. (13). 
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The perturbed current is given by the following expression: 


(a;/3Wk,a + V-k^aXjB) g (£'k) V/3X- 
k 

Let us now study the px coefficient, it is defined by the expression: 

Map = -^Tr'^{xpvi^^a + V]^,aXp)g{S]i) = 5(^)Tr [(*/3^k,a - XaVi^^p)S{g - H]^)]dg. 

We follow Smrcka and Streda approach and introduce two functions: 

dG^ dG^ 

Aap{g) = *Tr Wk,a-^Wk./3<5(7? - -ffk) - Vk,a5{r] - H]^)Vk,p-^ , 

Bapig) = *Tr [vk,aG%k,p5{g - H^) - fk.a<5(?? - Hu)vk,pG^] , 

where = [r] — ± i0+). Expression 






:Tr [a;a (G^ - G^) a;;? - XaXp (G^ - G^)] - (a o /?) 


is useful for further calculations. We also use the following expressions: 


• |a;a l^(G^)^ - (G^)^j Wk./jj - (a <-> ^) = 27riTr {XaVk^p - xpv^^a)-^Hv - Bk) . 


We finally obtain: 


^apiv) - \ ~ ^/3 “ ^«^/3 - (a ^ /5) (40) 


+ ^Tr {xaVk,p - xpVk,a) ^<5 iv - Hk) , (41) 

where the first term after integration over r] will result in a commutator [x^^xp], hence will vanish. Vanishing of 
this commutation relation will only happen when the whole Hilbert space is considered. The following identity for 
bounded spectrum is of use 

/ oo / -I jD ^(r^\\ roo r 

dd ^ J [^Kc^v^^pSiri - iJk) - - iJk)^k,^^J (42) 

= -i [ dg'XS[g- e„k] I Vk,a- -~7^^k./3 \ - (a /?) (43) 

n I [Snk-^^k] J„„ 

/ oo 

dgS [g - {Hk)nn\ (A,/3)„„j ("4k,a)„„ -{a^ 13) (44) 

-oo 

= z y] [(a^T^) {daTk)] - (a o /?) = 0. (45) 

n 

In the expression above (k) = i (^/ 3 Tk)j — {a ^ is the /c—space Berry curvature of the nth band. 

We observe a sum rule, ~ ^5 which will be uesful in further derivations. Another identity 


Badiv) = ^J2\ (dpT^)] 6 [t] - 5nk] - [a ^P) 
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is of use in further derivations. After all these transformations, one can show that 

/ rOO rOO rO r — OO \ / 

Map = dr) J dfi + J dr) J dfjj g{ri) {^Aap{v) - 


1 

2 dfj 


= dfi (^Aapiv) 


1 dBap{fl) 

2 dfj 


f dr]g{r]). 
Jo 


We then get a Berry curvature contribution to the response written as 

/ oo r / \ 1 

dfjS [fj - £„k] (SaTk) (dpTk) / g{v)dg - ggig) 

-oo LV / inn Jq 


— {a ^ P) 




kn 


(47) 

(48) 


(49) 

(50) 


where ci{x) = r]dr] — rj^^^^dr], in which a sum rule derived above in expression (42), namely (k) = 

0, can be used to eliminate the constant terms of the integral. Finally, the full expression for the response is given by 
two compact terms 


Ja- V 


(k)ci ^ ^ (^g^nk) (^;g^nk)^nk 9 (^nk) 


kn 


kn 




where we restored /cb factor and volume. 


(51) 


C. Torque response to temperature gradient 

In order to calculate the torque response to the temperature gradient, we define two contributions: 

(d^H) = ^ {dm [rpH + Hrp])^^ Vpx- (52) 

The first term is described by the Kubo formula, (^m^)ne ~ fhe second term is defined as = 

^ (^m we then formally rewrite the expression for torque as 

{dmH) = {Sm/3 + Mm/s) V^x = F^^V^X- (53) 

Calculations for the torque are similar to the ones presented for the particle current. The final result is given by the 
expression: 


(dmH) = - 


ksT J2 f^m^(k)ci (£nk) + J2^dmen k) (^;0^nk)^nk 2p 9 (^nk) 


kn 


kn 


V/3X, 


(54) 


where now f2[^^(k) = i |^(9mTk) (^/ 3 Tk)j — (m pp (5) is the mixed space Berry curvature of the nth band, and we 

restored the /cb and volume factors. 
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